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Abstract We study quantum spin systems in the ID, 2D square and 3D cubic 
lattices with nearest-neighbour XY exchange. We use the coupled-cluster method 
(CCM) to calculate the ground-state energy, the T = sublattice magnetisation and 
the excited state energies, all as functions of the anisotropy parameter 7. We consider 
S = 1/2 in detail and give some results for higher S. In ID these results are compared 
with the exact S = 1/2 results and in 2D with Monte-Carlo and series expansions. 
We obtain critical points close to the expected value 7 = and our extrapolated 
LSUBn results for the ground-state energy are well converged for all 7 except very 
close to the critical point. 
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1. Introduction and CCM formalism 

In this paper we consider the T = properties of the quantum spin system known 
as the XY-model, described by the Hamiltonian 

H = \j2K 1 + 7Ks? +P + (1 - 7K*W m the regime < 7 < 1 (1) 

i,p 

where index / runs over all N lattice sites with periodic boundary conditions, and p 
over the z nearest-neighbour sites. 

For s = 1/2 and ID this model was solved exactly by Lieb, Schultz and Mattis 
(1961) and its properties have been studied by many authors (see Niemeyer, 1967, and 
Barouch et al, 1971, for example). For higher spin in ID or in 2D (square) and 3D 
(simple cubic) useful results have been obtained using spin- wave theory (Zheng et al, 
1991), Monte-Carlo methods (Ding, 1992, Zhang and Runge, 1992), series expansions 
(Hamer et al, 1991) and, for 7 = 0, finite size extrapolations (Betts et al, 1996). 

In a recent paper (Bishop, Farnell and Parkinson, 1996, referred to as I), the 
coupled-cluster method (CCM) was applied to the XXZ model in the |A| < 1 regime. 
It was found that good results could be obtained by using a planar model state in 
which the spins are aligned in the xy— plane, as in the classical ground state, rather 
than along the z— axis. Here we shall use a similar model state for (1), again motivated 
by the classical ground state. 

For a description of the CCM applied to spin systems see Bishop et al. (1991) and 
also the references given in I. To calculate the ground state wave function \^/) of a 
spin system we start with a model state |<3>) and a correlation operator S such that 

= e 5 |$) 

For the Hamiltonian (1) we expect that in the ground state the spins are aligned 
in the xy— plane. We choose |<3>) to be a Neel state with spins aligned parallel and 
antiparallel to the x— axis. In ID this has the form 

|*> = |... < > < > < > < > < ► ...)• 

It is useful to introduce local axes such that each spin in |$) is pointing in the 
negative z— direction, by means of the following transformation: 

s x — > — s z , s y — > s y , s z — > s x left-pointing spins 
s x — > s z , s y — * s y , s z — > — s x right-pointing spins. 
Thus (1) becomes (with s ± = s x ± is y ) 



H = 9 Y,i As i s i+p + B (44 +P + sr^r+p) + c(^«r +P + sfs,y ] (2) 



2,P 



with 



A=-(l + 7 ), B 



-(l- 7 ), C=-5. 



For the correlation operator S 1 we choose a linear combination of creation operators 
relative to |<&), a creation operator being any combination of spin raising operators 
(s + in the local axes). Because of the form of (2) the total number of spin flips in 
each creation operator must be even. 

The simplest possible choice for S is to flip two spins, known as the SUB2 approx- 
imation scheme: 



where r runs over all distinct lattice vectors (r ^ for s = 1/2). 

The full SUB4 scheme involves 4-fl.ip configurations as well 2 flips and is too com- 
plicated to handle in general. However the most important extra term is the one with 
4 flips on adjacent sites. Including this term gives the SUB2+LSUB4 scheme which 
we have applied only in ID: 



A third approximation scheme is to include in S all possible combinations of spin 
flips within a region of size n, known as the LSUBn scheme. This is particularly 
useful for numerical extrapolation as a function of n, and will be discussed in detail 
in section 5. 

From the Schrodinger equation H\^f) = E\^f) we obtain the equation for the 
ground state energy: 



This equation is exact whatever approximations are made for S. 

To determine the coefficients b r and g± in the SUB2+LSUB4 scheme we operate 
on the Schrodinger equation with exp(— S) then one of the destruction operators and 
then by (<3>|: 




(3) 







($| S rw- s tfe 5 |$> = E[ 5 E^-r' + p - (A + ABh)b r 

V r' 



+ 2C6 r „ p + {B{2b\ + 2^ + 1) + A&!)5 Pir + Bg 4 S 3pjr ] = (6) 

($|s^s m s r+2 s r+3 e ~ 5 # e5 | $ ) = A(6 2 + 36^ + 26x63) - 45(6x6264 + hb 2 3 + b 2 2 b 3 ) 

- 4C(26x6 2 + 6 2 6 3 ) + g A [B(2b 5 - 26 3 - 860 - A] = (7) 

The corresponding equations for the SUB2 scheme are obtained by setting g± = 
everywhere in the first of these and ignoring the second. 

These coupled non-linear equations are solved by first Fourier transforming Eq. (6) 
and then solving the resulting equations and Eq.(10) self-consistently. For dimension 
d we obtain: 

T(q) = £ e ir *b r , b r = r d d q(27r)- d e-^T(q), 7(a) = - £ <*« 

r J —7V Z p 

61 = [ W rf d g(27r)- d 7 (g)r(g), X ± = £ 6 r 6 r+p = f d d g(27r)- d 7 (g)r 2 (g) 
leading to 

aT 2 {q) + 6r(g) + c = 0, 

where 

a = S 7 (g), 6= -A-4fi6 1 + 2C 7 (g), 
c = [5(26? + 2^4 + 1) + A6i)]7(g) + %4 7 (3g) - BXi - 2C6x 
with the usual solution 

, -6 + V6 2 - 4ac 
r(?) = 2^ • 

The equations can now be solved numerically in a self-consistent way. 

Results for the ground state energy using the SUB2 and SUB2+LSUB4 approx- 
imation schemes are shown in Figures 1 and 3, and Tables 1 and 2 for ID and 2D. 
The LSUBn results are discussed in section 5. 

A notable feature of the CCM is the existence of terminating points as a function 
of 7. These are believed to correspond to the actual T = phase changes, known 
to be at 7 = in ID and believed also to be at 7 = on symmetry grounds for 
2D and 3D. In ID terminating points only occur if correlations of infinite range are 
explicitly included in S and occur at 7 = —0.10789 in SUB2 and at 7 = —0.09605 in 
SUB2+LSUB4. In 2D there is a terminating point at 7 = -0.03033 in SUB2. These 
are reasonably close to 7 = considering the simple nature of these approximations. 



In 2D and 3D terminating points can also occur within the LSUBn scheme as 
described in section 5. 

2. In-plane Sublattice Magnetisation 

In the CCM the bra ground state is not in general the Hermitian conjugate of the 
ket state. Instead we introduce a new operator S such that 

(*| = (Q\Sexp(-S) . 
The SUB2+LSUB4 approximation for S is 

S = 1 + Yl (l, S ~ b r S i S l+r + hs^S^S^S^ (8) 

where r runs over all distinct lattice vectors (with r 7^ for s = 1/2). 

The bra-state equations are found variationally by taking the partial derivatives 

of 

H = (V\H\V) 

with respect to the ket-state coefficients. By CCM theory (Bishop et al. 1991) these 
derivatives must be equal to 0. Hence we obtain two bra state equations: 



~K7~ = NY\2BY,~br>br>-r+ P - (A + ABh)b r + 2Cb r - p 
+(B + (A + 4S6i)6i - 45 ^ & r /&r')<W 

r> 

+g 4 /2{[AA(b 1 + 63) - 85(6 2 6 4 + b\) - lQBg 4 - 16C6 2 ]5 p , r 

+ [12Ab 2 - 85(6x64 + 26 2 6 3 ) - 867(26! + b 3 )]5 2p , r 

+ [4Ah - 85(26x6,3 + b\) - 4Bg 4 - 8Cb 2 ]5 3p , r 

-856x6 2 5 4p ,r + 45(?4<W}] = (9) 



dH 



N 



5(26x + 63) + ^(25(65 - 63 - 46i) - A) 







(10) 



Again we perform a Fourier transform on Eq.(9) and the resulting equations and 
Eq.(10) may be solved self-consistently in order to obtain the bra-state correlation 
coefficients. 



Finally the results are used to calculate the magnetisation using the formula for 
SUB2: 



M = -2(V\sf\V) 

and for SUB2+LSUB4: 



= 1 - 2 ^ KK 

r 



M = -2<*|sf = 1 - 2]TMv - 8^4 



3. Excitations 

A similar method can be used for the excited state energies, introducing the oper- 
ator 

Xi = XjSj , i belongs to one sublattice only 

i 

leading to 



(<!>\sre- s [H,X l }-e s \^} = ~z(A + ABb l )X l + B £ b r X l+r+p = e x X x (11) 

r,p 

and hence, via Fourier transform 

=> e(q) = -\z{A + ABb x ) + 5z 7 (g)r(g) (12) 



4. General spin s, (SUB2 only) 

We have also considered the general case of s > 1/2 within the SUB2 approxima- 
tion scheme. The main features are as follows. 

The correlation operators S and S are the same as before (without g±). The ket 
state equations become: 



j r —p 



($\s-[sr +r e- s He s \<5>) = As 2 ^ [As 2 B £ b r ,b r - r , +p - 2s(A + 4B6i)6 r + AsCb r 

P r' 

+{B{2b\ + 1) + Ab 1 )5 p>r ] = (13) 

and the energy is 

($\e- s He s \<5>) = 2s 2 zN (^A + ^B^j (14) 



Using these equations we find for the ground-state energy per spin of the s = 1 
system at the 7 = point the value —1.09179. This compares with a numerical result 
from extrapolating rings with N < 14 of —1.1157 ± 0.0003. This is a very similar 
accuracy to that obtained using SUB2 for s = ~ at the same point. 

There are similar modifications to the bra state equations which become 



OH 

— = 4s 2 Nj2[8s 2 Bj2br'b r '-r+p-2s{A + ABb 1 )b r + 4sC6 r _ p 

+ (B + (A + 4B&i)6i - 8sBj2br'b r ')S P) r] =0 (15) 

r' 

Finally the magnetisation is given by: 

M= — (*|sf|*) = 1 -4s^M> r (16) 



s 



r 



5. The LSUBn Approximation 

The LSUBn scheme contains all possible (connected and disconnected) terms in 
S which are contained within a 'locale' of size n. We use all possible connected 
configurations of n spins to define this locale; in ID we may see that this locale 
is simply a chain of length n. Disconnected and connected configurations of less 
than n spins are then generated by successively removing sites from the original 
connected configurations of n spins, thus covering all possibilities. The lowest order 
LSUBn approximation scheme is the LSUB2 (i.e., SUB2-2) approximation in which 
only a single nearest-neighbour, two-body term is retained in S. We note that the 
Hamiltonian of Eq. (2) includes products of the spin operators which contain even 
numbers of these spin operators only. This means that the ground state contains only 
even numbers of spin flips with relation to the model state. We restrict the LSUBn 
approximation to include only those configurations which contain an even number of 
spin raising operators. A further restriction is that each fundamental configuration 
must be independent of all others under the symmetries of both the lattice and the 
Hamiltonian; we note that both lattice and the Hamiltonian have identical symmetries 
for the XY model. 

Tables [I], 0, and [$] show the numbers of fundamental configurations for given 
LSUBn approximation level, and we can see from these tables that the number of 
configurations grows very rapidly with n. Hence, for higher-order approximations we 
need to enumerate all possible configurations computationally, and we furthermore 
need to obtain and solve the CCM LSUBn equations computationally also. A full 
explanation of the computational method used here is given in Zeng et al. (1997). It is 
now possible to obtain values for the ground-state energy and sublattice magnetisation 
for the LSUBn approximation scheme. Results for these quantities are given in Figs. 



1,2,3 and 4, and results at the isotropic point of 7=0 are given in Tables [I], [| and 
|3|. A simple extrapolation of the ground-state energy and sublattice magnetisation 
has also been carried out by plotting the ground-state energy against 1/n 2 and the 
sublattice magnetisation against 1/n, and then performing polynomial fits on this 
data. The extrapolated LSUBoo results obtained from this simple, 'naive' approach 
are shown in Tables [1], 0, and [3|. The results are clearly at least as good as obtained by 
series expansion. Results in 2D and 3D are especially valuable since no exact results 
are available. 

Another consequence of this approximation scheme is that the second derivative of 
the ground-state energy is found to diverge for some critical value of the anisotropy 
parameter, denoted 7 c (n), in 2D and 3D only. These points are related to phase 
transitions of the true ground state of the system (Zeng et ai, 1997), and the results 
for given LSUBn approximation level are shown in Tables |2|, and |3|. We note that 
critical j c (n) approaches 7=0, the point at which the true phase transition point 
is believed to be (in all dimensions), with increasing approximation level. Again, 
a simple extrapolation of the LSUBn critical points is attempted by plotting 7 c (n) 
against 1/n 2 , as in Bishop et al. (1994), and the extrapolated results are also shown 
in Tables ^| and |3|. 
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Figure captions 

Fig.l 

Results for the CCM ground-state energy of the one dimensional XY model. The 
terminating points of SUB2 and SUB2+LSUB4 schemes are indicated. 

Fig-2 

Results for the CCM ground-state sublattice magnetisation of the one dimensional 
XY model. 

Fig-3 

Results for the CCM ground-state energy of the square lattice XY model. All the 
approximation schemes have terminating points except LSUB2. 

Fig.4 

Results for the CCM ground-state sublattice magnetisation of the square lattice XY 
model. 



Table 1: Ground-state energy and sublattice magnetisation for the one dimensional 
XY model at 7 = compared to exact results of McCoy (1968). Nf indicates the 
number of fundamental configurations for a given LSUBn approximation level. 



LSUBn 


N/ 


E 9 /N 


M 


LSUB2 


1 


-0.303813 


0.837286 


SUB2 




-0.310377 


0.779517 


LSUB4 


4 


-0.314083 


0.722916 


LSUB6 


13 


-0.316301 


0.660064 


LSUB8 


43 


-0.317137 


0.617624 


LSUB10 


151 


-0.317542 


0.586067 


LSUBoo 




-0.31829 




Exact 




-0.318310 


0.0 



Table 2: Ground-state energy and sublattice magnetisation for the square lattice XY 
model at 7 = compared to series expansion calculations of Hamer, Oitmaa and 
Zheng [2]. Nj indicates the number of fundamental configurations for a given LSUBn 
approximation level, and also shown are the critical values of 7 for the anisotropic 
model - where the value in parentheses is the estimated error in the final decimal 
place shown. 



LSUBn 


N/ 


E s /N 


M 


7c M 


LSUB2 


1 


-0.540312 


0.949634 




SUB2 




-0.546325 


0.918953 


-0.030(1) 


LSUB4 


10 


-0.547267 


0.915768 


-0.175(1) 


LSUB6 


131 


-0.548329 


0.901357 


-0.073(1) 


LSUB8 


2793 


-0.548616 


0.893665 


-0.04(1) 


LSUBoo 




-0.54892 


0.869 


0.00(1) 


Series Expansion 




-0.5488 


0.872 


0.0 



Table 3: Ground-state energy and sublattice magnetisation for the cubic lattice XY 
model at 7 = 0. N/ indicates the number of fundamental configurations for a 
given LSUBn approximation level, and also shown are the critical values of 7 for 
the anisotropic model - where the value in parentheses is the estimated error in the 
final decimal place shown. 



LSUBn 


N/ 


E g /N 


M 


7c M 


LSUB2 


1 


-0.786866 


0.971488 




SUB2 




-0.790901 


0.958282 


-0.01666(1) 


LSUB4 


13 


-0.791224 


0.958648 


-0.172(1) 


LSUB6 


327 


-0.791702 


0.954759 


-0.071(1) 


LSUBoo 




-0.79201 


0.948 


0.01(1) 



E g /N 



-0.53 



-0.58 



-0.63 




-0.20 



0.20 



1.00 



0.95 



0.90 



0.85 



0.80 



0.75 




0.70 



-0.5 



0.0 



-- LSUB2 
-- LSUB4 

— LSUB6 

— SUB2 

— LSUB8 

0.5 



1.0 



Y 



